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There is recent evidence that some SiC X grains from the Murchison meteorite may contain 
supernova-produced ^-process n B and or 7 Li encapsulated in the grains. The synthesis of X1 B 
and 7 Li via neutrino-induced nucleon emission (the (/-process) in supernovae is sensitive to the 
neutrino mass hierarchy for finite sin 2 2(9i3 > 0.001. This sensitivity arises because, when there is 13 
mixing, the average electron neutrino energy for charged-current neutrino reactions is larger for a 
normal mass hierarchy than for an inverted hierarchy. Recent constraints on #13 from the Daya Bay, 
Double Chooz, MINOS, RENO and T2K collaborations all suggest that indeed sin 2 20 13 > 0.001. 
We examine the possible implications of these new results based upon a Bayesian analysis of the 
uncertainties in the measured meteoritic material and the associated supernova nucleosynthesis 
models. We show that although the uncertainties are large, they hint at a marginal preference for 
an inverted neutrino mass hierarchy. We discuss the possibility that an analysis of more X grains 
enriched in Li and B along with a better understanding of the relevant stellar nuclear and neutrino 
reactions could eventually reveal the neutrino mass hierarchy. 



PACS numbers: 14.60.Pq, 26.30.-k, 25.30. Pt, 97.60.Bw 

I. INTRODUCTION 

Oscillations in the three-neutrino flavor mixing sce- 
nario are described by three angles 812, #23, and #13 plus 
a CP-violating phase 6c p- At the present time solar, 
atmospheric, and reactor neutrino oscillation measure- 
ments [THS] have provided information on the neutrino 
mass differences, i.e Am 2 2 = \m\ — m 2 | — 0.000079 
eV [5] and Am 2 3 w |Am 2 3 | « 2.4 x 10" 3 eV 2 [1 ©. 
These measurements, however, are unable to determine 
the mass hierarchy, i.e. whether Am 2 3 > (normal) 
or Am 2 ^ < (inverted) is the correct order. Also, so- 
lar, atmospheric and reactor neutrino oscillation experi- 
ments [IHS] have determined 6*12 and #23 to reasonable 
precision. However, only recently have measurements 
of #13 become available. The three best best current 
measurements [with 8 c p = and 6*23 = 7r/4] are from 
the the Daya Bay experiment that has recently reported 
[7] sin 2 26»i3 = 0.092 ± 0.016(stat) ± 0.005(syst), the in- 
dependent recent report [8] from the RENO collabora- 
tion of sin 2 20i3 = 0.113 ± 0.013(stat.) ± 0.019(syst.) and 
the value [S] from the Double Chooz collaboration of 
sin 2 26<i3 = 0.086 ± 0.041(stat.) ± 0.030(syst.). 

These results are mutually consistent and also consis- 
tent with the previously reported upper limit from the 
MINOS collaboration [B] of sin 2 26»i 3 < 0.12(0.20) for the 
normal (inverted) hierarchy, and with the 90% C.L. lower 
limit and upper limits to #13 from the T2K collaboration 
[I0] of 0.03(0.04) < sin 2 26*13 < 0.28(0.34). As encour- 
aging as these new results are, however, the data do not 



yet determine whether the normal or inverted hierarchy 
is the correct ordering of mass eigenstates. 

In this context we note that previous studies [TTUT5] of 
the nucleosynthesis of the light isotopes via the (/-process 
in supernovae have pointed out that for a finite mixing 
angle 6*13 > 0.001 the relative synthesis of 7 Li and n B 
in the ^-process is sensitive to the mass hierarchy. Since 
6*13 > 0.001 is indeed implied by the Daya Bay +RENO 
+ Double Chooz results to better than the 5a C.L., it 
is worthwhile to reconsider the supernova (/-process as a 
means to constrain the neutrino mass hierarchy. More- 
over, there has also recently appeared a possible discov- 
ery |16) of (/-process supernova material in SiC X grains 
from the Murchison meteorite. The goal of this paper 
is to examine whether constraints on the neutrino mass 
hierarchy can be determined from these data. 

It has recently been pointed out |17j . however, that the 
supernova models themselves can lead to variations in the 
synthesis of 7 Li and n B that are much larger than the 
anticipated mass hierarchy effect due to uncertainties in 
the stellar thermonuclear reaction rates. Moreover, the 
previous studies only considered a single 16.2 M super- 
nova model, while there is a significant dependence of the 
(/-process yields on supernova progenitor mass which can 
vary 

[H1EI] from 10 to 25 M Q . Nevertheless, rather than 
to give up on this quest, one should keep in mind that 
there are prior constraints on the degree to which the nu- 
clear reaction rates can be varied, and on the probability 
for any particular supernova progenitor model. Our goal 
here is to consider the likelihood of one neutrino mass 
hierarchy over another in a statistical analysis that takes 
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proper account of all of these effects. 

Indeed, the formulation of Bayesian statistics [20] [2T] 
provides just such a framework in which to analyze the 
probability that any particular model is true even in the 
context of uncertain data and large variations in under- 
lying predictions due to model parameter uncertainties. 
Moreover, this approach provides a means in which to in- 
corporate all prior constraints on model parameters and 
input data. In the remainder of this paper, therefore, we 
review the underlying model predictions, and their un- 
certainties along with an analysis of the meteoritic con- 
straints. These are then applied in a five dimensional 
Bayesian analysis. We find that in spite of the large un- 
certainties, there remains a marginal preference for an 
inverted hierarchy at the level of 74%/26% compared to 
a prior expectation of a 50%/50%. 



II. THE ^-PROCESS 

The ^-process [551 122] is believed to occur in core- 
collapse supernovae. As the core collapses to high tem- 
perature and density, neutrinos of all flavors are gener- 
ated thermally in the proto-neutron star formed by the 
collapse. These neutrinos emerge from the proto-neutron 
star at relatively high temperature such that some neutri- 
nos can participate in neutrino-induced nucleon emission 
reactions. The effects of these neutrino interactions are 
usually negligible. However, the i^-process can be a ma- 
jor contributor to the production of very rare isotopes 
not produced by other means such as the light isotopes, 
7 Li, 11 B, 19 F, and the heavy rare odd-odd isotopes, 138 La 
and 180 Ta [TH1I13H1S]. 

The isotopes 7 Li and n B, are of particular interest 
for the present work. These nuclides can be copiously 
produced (e.g. ( 7 Li/H)/( 7 Li/H) Q ~ 10 4 [Mil]) in su- 
pernovae as the neutrinos flow through the outer He layer 
of the exploding star. Although 7 Li can be also be pro- 
duced in the Si shell [23], it will be quickly destroyed by 
photodisintcgration. Therefore, in the present work we 
expect the final ejected abundances to be predominantly 
from the He shell and can disregard production in the Si 
shell. Key to the present context is the sensitivity [15] 
of the produced 7 Li/ 11 B abundance ratio to the neutrino 
mass hierarchy. This arises because neutrino oscillations 
can increase the temperature characterizing the spectrum 
of reacting neutrinos as they transport through the SN 
ejecta [24] , This can increase the rates of (/-reactions 
relative to models without oscillations thus affecting the 
yields of the light elements. 

The dependence on #13 enters because there is a res- 
onance of the 13-mixing in the outer C/O layer [TTMl5| . 
The associated increase in the electron neutrino temper- 
ature increases the rates of charged-current z/-process re- 
actions in the outer He-rich layer. The heavy i/-process 
elements are produced in the O-rich layers below the res- 
onance region [TS] and are thus not affected by the 13 
neutrino oscillations. The light isotopes 7 Li and n B, 



however, are significantly produced in the outer He layer 
as described below, and thus are sensitive to the 13 mix- 
ing. Moreover, the produced 7 Li / 11 B ratio depends not 
only on the mixing parameter, #13, but also on the neu- 
trino mass hierarchy if the mixing angle is as large as 
suggested by current measurements. 

The detailed dependence on the neutrino mass hier- 
archy is as follows [TTI - H5] . For a normal hierarchy, 
the adiabatic 13-mixing resonance for neutrinos in the 
O/C layer causes the v e energy spectrum in the He/C 
layer to be comparable to that of the in the O-rich 
layer. Without oscillations, 7 Be is produced through 
the 4 He(z/, z/n) 3 He(a, 7) 7 Be reaction sequence. How- 
ever, when there are neutrino oscillations, the rate of 
the 4 He(V e , e~p) 3 He reaction exceeds the 4 He(V, j/n) 3 He 
rate so that the production of 7 Be in the He layer is much 
larger if neutrino oscillations occur. 

The production of 7 Li is also greater when neutrino os- 
cillations occur. In this case, however, the enhancement 
is much less than that for 7 Be. The main production 
path for 7 Li is the 4 He(z/, v'p) 3 H(a, 7) 7 Li reaction se- 
quence and the corresponding charged-current reaction 
is the 4 He(z/ e , e + n) 3 H reaction. There is, however, no 
resonance for the antineutrinos, so this rate is insensitive 
to the 13 mixing. 

Neutrino oscillations also affect the production of 
n B and U C similarly to that of 7 Li and 7 Be. Dur- 
ing the i/-process, n C is normally produced via the 
12 C(V, z/n) n C reaction, plus a smaller contribution from 
the 12 C(v e ,e~p) 11 C reaction. When oscillations occur, 
the rate of the 12 C{y e , e~p) 11 C becomes larger than the 
rate without oscillations by about an order of magnitude 
so that the production of 11 C in the He layer is larger 
when neutrino oscillations occur. 

The abundance of n B with oscillations, how- 
ever, is only slightly larger when oscillations oc- 
cur. The main production process for U B is 
the 4 He(t/, ^'p) 3 H(a, 7) 7 Li(a, 7) 41 B reaction sequence. 
The corresponding charged-current reaction is the 
4 He(i? e ,e + n) 3 H reaction. About 12%-16% of the n B 
in the He layer is also produced from 12 C through the 
12 C(v, z/p) n B and the 12 C(D e , e+p) n B reactions. How- 
ever, oscillations produce no increase in the n B produc- 
tion through the 12 C(P e , e + n) n B or 4 He(f'e, e + n) 3 H re- 
actions, because the antineutrinos have no resonance. 

In the case of an inverted hierarchy, the production of 
7 Li and n B is greater than for a normal hierarchy. The 
4 He(D e , e + rt) 3 H and 12 C(D e , e + n) n B reaction rates are 
larger because of an adiabatic resonance of v e -f-> v^ T . 
However, the increased production of 7 Li and n B is less 
pronounced than in the case of a normal hierarchy. This 
is because the average v e energy is already greater than 
the v e energy at the neutrino sphere. Therefore, the av- 
erage energy is not changed much by the oscillation with 
u llT (D llT ). The mass fractions of 7 Be and 11 C are slightly 
larger in the outer parts of the star, and slightly smaller 
inside the He layer. There is, however, no 13-mixing res- 
onance for v e . Hence, there is no significant conversion 
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of v e -f-s- v^ T . At the same time, some of the produced 
7 Be and 11 C is affected via capture neutrons from the 
4 He(z/ e , e + n) 3 H reaction. 

The final result from all of the above is that the emer- 
gent 7 Li I/ / 11 B I/ ratio from the He layer is sensitive to 
whether the neutrino mass hierarchy is normal or in- 
verted [I3TH5] . (We hereafter denote isotopes produced 
in the ^-process by Z„.) 

Figure [T] shows an extended plot of the variation of 
the 7 Li y / By ratio as determined in [15] as a function of 
sin 2 (26»i 3 ). The key point is that for sin 2 (26»i 3 ) > 1CT 3 
there is a significant difference in this ratio between the 
normal and inverted hierarchies due to the effects of a 
higher v e energy in the normal hierarchy. 




FIG. 1: Produced 7 Li/ J1 B abundance fTS] as a function of 
mixing angle for both a normal and inverted neutrino mass 
hierarchy. This ratio varies in models with different neutrino 
temperatures in the range indicated by the lower and upper 
solid lines. The width of the shaded region indicates the 2a 
confidence limits to sin 2 (2^13) from the Daya Bay result [7|. 
The top of the shaded region is the 2a upper limit on the 
observed z/-process 7 Li/ 11 B ratio as deduced here from the 
SiC X grains. 

The results of Refs. [T3HH], however are only based 
upon a single 16.2 M model and did not account for 
variation of ^-process 7 Li and n B yields with progenitor 
mass [IB] and uncertainties in underlying thermonuclear 
reaction rates. Indeed, it has been demonstrated [T7] 
that there are large uncertainties in the boron produc- 
tion due to the fact that uncertainties in the 3a -^ 12 C 
and 12 C(a, 7) 16 nuclear reaction rates can lead to dif- 
ferent core configurations and ultimately, differences in 
the produced U B and 7 Li by factors larger than the sen- 
sitivity to the neutrino mass hierarchy. 

What is needed, therefore, is a proper statistical anal- 
ysis of the uncertainties in all model parameters in or- 
der to determine whether ^-process nucleosynthesis is a 



possible means to constrain the neutrino mass hierarchy. 
The purpose of this paper is to explore such a possibility 
based upon a Bayesian analysis of the underlying uncer- 
tainties in the supernova models and measured meteoritic 
material. 



III. METEORITIC //-PROCESS MATERIAL 

All of this is well and good, but would be of no value 
without a measurement of ^-process 7 Li and n B in super- 
nova material enriched by the ^-process. The detection of 
such rare isotopes directly in supernova remnants would 
be exceedingly difficult. Fortunately, however, there is 
another way. For many years it has been recognized that 
SiC grains in carbonaceous chondrite meteorites can con- 
tain a component of material freshly synthesized in a 
supernova and then trapped in the grain. Most of the 
pre-solar grains appear to represent the ejecta from AGB 
stars of 1-3 M ; however a small fraction (~ 1%) of them 
(the so-called X grains) can be attributed to ejecta from 
core collapse supernovae [26 28 . The isotopic ratios of 
these X grains exhibit 12 C/ 13 C higher than solar and 
14 N/ 15 N lower than normal. They also have enhanced 
28 Si. These isotopic features are all characteristic of core 
collapse supernovae. Moreover, they contain decay prod- 
ucts of radioactive isotopes synthesized in supernovae 
such as 26 Al (t 1/2 = 7 x 10 5 yr) and 44 Ti (t 1/2 = 60 
yr), establishing that these grains indeed formed from 
encapsulated fresh supernova material [28] . 

In this context it is of particular note that a recent 
study [16] of 1000 SiC grains from a 30 g sample of the 
Murchison CM2 chondrite gS] found 12 X grains that 
show resolvable anomalies in Li and/or B. In particular, 
an average of the 7 best-case single grains indicated a 
2a detection of excess n B relative to 10 B (i.e. n B/ 10 B 
= 4.68±0.31 compared to the solar ratio of 4.03 [501131]) 
along with an enrichment of Li/Si and B/Si relative to 
solar. As 11 B is expected to be the main product of light 
element ^-process nucleosynthesis [TS], this hints at an 
average enrichment of ~ 15% of ^-process n B. On the 
other hand they determined a slightly below solar aver- 
age ratio of 7 Li/ 6 Li = 11. 83 ±0.29 compared to the mete- 
oritic solar ratio of 12.06 [32] , However, the uncertainty 
is consistent with some 7 Li excess even at the la level. 
Hence, this result can be used to set upper limits on the 
7 Li„/ 11 B„ production associated with these grains. 

Assuming that the observed grain 7 Li and n B abun- 
dances can be decomposed into average solar plus re- 
process contributions, we then write: 
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and similarly, 
io B 



11 



4.03 



4.68 ±0.31 



(2) 



where we explicitly note that neither b Li nor U B is 
significantly produced in the ^-process |15j . Next we 
adopt the solar ratios 7 Li/ 6 Li = 12.06 ± 0.02 [32] and 
n B/ 10 B= 4.03 ± 0.04 [30] consistent with carbonaceous 
chondrites. We note that other 7 Li/ 6 Li ratio values ex- 
ist in the literature based upon a combination of solar 
photospheric and meteoritic data, e.g. 7 Li/ 6 Li= 12.176 
[S], 12.33 [33], 12.177 [34]. However, the value adopted 
here is the one found in CI meteorites and hence the best 
representation of primitive solar system material. 
From Eqs. (1) and (2) we can deduce 



r Li u 



Li 



-0.35 ±0. 



/ x 



(3) 



where the subscript X denotes the 6 Li/ 10 B ratio appro- 
priate to the X-grain samples under study. Although 
it is probably safe to adopt the solar isotopic ratios for 
the grain contaminant, elemental ratios are not appro- 
priate as there is no guarantee that the non-r/-process 
Li/B in the X grains is identical to the solar value due 
to the different volatilities of Li and B. Fortunately, how- 
ever this ratio can be determined directly for the grains 
from the measured elemental ratios in |16j . The average 
for the 7 single grains is Li/Si= 8.27 ± 0.06 x 10~ 5 and 
B/Si= 4.13 ±0.11 x 10~ 5 . 
From this we write: 



fLi 0+ 7 Li,)/s L i + l - 2 00 ±0 04^ 
("B +» B„)/WB + l\ x - 2m±0m { 6 Li G 

"(4) 

Then inserting the X grain measured isotopic ratios 
given in Eqs. (1) and (2), we deduce an average X-grain 
ratio of ( 6 Li / lo B Q ) x = 0.89 ± 0.06. Inserting this into 
Eq. (3) leads to 7 LL,/ n B„ = -0.31 ±0.42 and the follow- 
ing upper limits on the X-grain ^-process isotopic ratio: 



n B„ 



< 0.53(2cr 95% C.L.) 

< 0.95(3cr 99.7% C.L.) 



(5) 



The width of the shaded region on Fig. [T] shows the 2a 
(95% C.L.) determination of sin 2 (26*13) from the Daya 
Bay analysis. The top of the shaded regions is the 2a 
upper limit to the 7 Li J/ / 11 B J/ ratio from the meteoritic 
SiC X grain analysis. 



IV. BAYESIAN ANALYSIS 

Even though the meteoritic data only provides an up- 
per limit, one can use Bayesian statistics to quantify the 
degree to which the information in this evidence supports 



an inverted or normal hierarchy neutrino mass hierar- 
chy. In particular, one can use Bayes' theorem [20] [21] 
to ascertain the probability P(M,i\D) that a given model 
Mi {i ^inverted or normal mass hierarchy) is true based 
upon a data set D (the meteoritic evidence). 
According to Bayes theorem we have: 



P(Mi\D) 



P{D\M t )P{Mj) 



(6) 



where P(D\Mi) is a likelihood function giving the proba- 
bility of the data to be reproduced by model M, (includ- 
ing all uncertainties in the model parameters) and P(Mi) 
is the prior probability (e.g. for our purposes we take to 
be 50/50) for an inverted versus normal hierarchy. The 
sum in the denominator is over all possible models, which 
in our case is just the inverted plus normal neutrino mass 
hierarchy. We will ignore here the possibility of no oscil- 
lations as sin 2 26*13 = is now ruled out at the level of 
5.2cr [7]. 

The key to this analysis is to identify the likeli- 
hood functions P{D\Mj), To do this we separate the 
likelihood into the probability distribution P(a k \Mi) of 
the parameters a k of the model and the probability 
P(E, Z, D\Mi, a k ) that a given model i with a set of pa- 
rameters a k produces a result Z that agrees with the ob- 
served data D that has an experimental error E. Then 
we have 

P(D\Mi) = J dEdZda k P{E,Z,D\M l ,a k )P{a k \M l ) , 

= J dEdZda k P(D\M i ,a k ,E,Z)P(Z,E\M i ,a k )P(a\M i ) .(7) 

Now for each set of parameters a k there is a unique 
prediction Zi for model Mi, so P(Z\Mi,a k ) — S[Z — 
Zi(Mi,a k )]. Moreover, the probability that a given Zf 
agrees with the observed data plus error D + E is given 
by the error distribution in the data (taken here to be 
gaussian). For the present application D ± erg = —.31 ± 
0.42 as deduced from the meteoritic data, so that we have 



P{E\Zi(Mi,a k )) = 



2naE 



exp[-(Z,-C) 2 /2a|) . (8) 



Next we need to specify the probability distribu- 
tions P(a k \Mi) for the parameters a k of the mod- 
els. For our purposes we identify 5 parameters, none 
of whose prior probability distributions depend upon 
whether the hierarchy is normal or inverted. So Eq. ^ 
becomes a five-dimensional integral over the probabil- 
ity distributions of each parameter weighted by the 
gaussian experimental error. These parameters are: 

Sin 26*13, T v , i?3a ; R\2Cai Mp rQ g as listed in Table |lj 

The parameter likelihood functions were decided in 
the following way. We adopt the Daya Bay result [7] 
of sin 2 26*13 = 0.092 as it has the smallest systematic 
error and is close to the weighted average of the three 
best recent determinations. Since the Daya Bay result is 
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dominated by statistical error, it is sufficient for our pur- 
pose to simply add the statistical and systematic errors 
in quadrature. We then adopt an overall gaussian proba- 
bility distribution for sin 2 2013 with a standard deviation 
of a = 0.017. 

For the rate of the 3a — > 12 C reaction, as in [T7] we 
consider an overall multiplier R 3a , where a value of unity 
corresponds to the standard evaluation [35] with a 12% 
gaussian error. Similarly, we consider a multiplier Rci2a 
for the 12 C(a,7) 16 reaction, with a best value [30] of 
Rci2a = 1-2 ± 0.25 times the evaluation in Ref. [55] . 

The probability distribution for a progenitor star to 
have a mass m is given by the initial mass function [37] 
cj)(m) = m~ x , with x = 2.65, and the normalization can- 
cels out in this and all parameter probability distribu- 
tions. The range of progenitor masses that can undergo 
core-collapse supernovae, however is limited to a range of 
10 to 25 M [IH]. Below 10 M the star does not form 
an iron core. Above 25 M a black hole probably forms 
rather than a neutron star. This interrupts the flux of 
neutrinos from the core so that there is no (/-process. 

For the neutrino temperature, the situation is a lit- 
tle more complicated. There are in fact six tempera- 
tures corresponding to the electron, muon, and tau neu- 
trinos plus their anti-particles. Although the cross sec- 
tions are roughly proportional to the neutrino temper- 
ature squared, this scaling cancels to first order when 
taking the ratio of yields. Even so, there is a small resid- 
ual sensitivity of the final 7 Li/ 11 B ratio to the neutrino 
temperatures which must be included. 

Fortunately, one can constrain phcnomenologically the 
relative neutrino temperatures and their possible vari- 
ations. For example, a recent detailed analysis [25] of 
the relative production of 180 Ta and 138 La in the v- 
process fixes the temperature of the electron neutrinos 
and their anti-particles to be around 4 MeV. On the other 
hand, successful r-process nucleosynthesis conditions re- 
quire a neutron-rich environment in the SN neutrino- 
driven winds. This requires that the electron neutrino 
temperature be slightly lower than that of the anti- 
electron neutrino temperature [TT], i.e. T Ve — 3.2 MeV 
for T Pe 4.0 MeV. Additionally, the observed galactic 
chemical evolution of boron and the meteoritic 11 B/ 10 B 
ratio have been shown 12J to constrain the and v T 
neutrino temperature to be about 5 MeV with an al- 
lowed range of T v = 4.3 to 6.5 MeV. This corresponds 
to Model ST in Ref. [T5]. We therefore adopt these as 
the best set of average neutrino temperatures. We then 
consider the range of models (1,2,LT and ST) studied in 
|15j as a suitable representation of the range of variation 
of the 7 Li/ 11 B ratio with neutrino temperature in core 
collapse models. These models incorporate the follow- 
ing ranges: T Ue = 3.2 - 4.0 MeV; T Pc = 4.0 - 5.0 MeV; 
T Ufl T = 5.0 - 6.5 MeV. The span of 7 Li/ u B ratios from 
this range is represented by the two sets of solid lines 
for each scenario in Figure [T] We treat this as a top-hat 
distribution in the integration of Eq. Q. The sensitiv- 
ity of the neutrino process to these models and also to 



TABLE I: Parameter likelihood functions P(a,k\Mi) 
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x = 0.92 
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T v (MeV) 
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T v = 3.2 - 6.5 


(see text) 


m 



variations in the neutrino cross sections has been care- 
fully studied in [15] and does not significantly impact the 
discussion here. 

Another sensitivity arises from the swapping of neu- 
trino energy spectra due to the self-interaction among 
neutrinos [38H4T] . It is well known that no one has 
yet succeeded in solving exactly the many-body prob- 
lem of three flavor oscillation due to the neutrino self- 
interaction. Nevertheless, as is demonstrated by Fogli et 
al. [ID], several cases of complete swapping of the high 
energy neutrinos between e- and fi- plus r-neutrinos (or 
between their anti-particles) should mimic the flavor os- 
cillation induced by an adiabatic MSW high-density res- 
onance: See Figs. 5 and 8 of Ref. [30] • This illustrates 
the maximum effect on the charged current interactions 
of the neutrino-process nucleosynthesis of interest to the 
present study. Hence, even should that effect be included 
in our adopted neutrino temperatures, the present dis- 
cussion would not change drastically. This is because the 
swapping effect is expected to result in the same isotopic 
ratios, i.e. 7 Li/ 11 B ~ 0.8 for a normal hierarchy and 
~ 0.6 for an inverted hierarchy even for a small mixing 
angle #13 < 0.001 which is excluded by current experi- 
ments. 

So now the likelihood for an inverted vs. a normal 
hierarchy becomes a five dimensional integral over the 
parameters listed in Table [I] 



P{Mi\D) 



J d 5 a k 



,-(D-Zif/2o 



zIL k P(a k \Mi 



-(D-Z z ) 2 /2a 



E U k P(a k \M t ) 



(9) 



The final task is then to specify the model prediction 
Zi(a,k, M^. In this we are aided by the fact that the 
16.2 M Q model calculation of [H] agrees well with the 15 
M Q models of [THJ [TO] and also those of 17 . Moreover, 
the relative effect of the oscillations on the 7 Li/ 11 B ratio 
depends only on the mixing angle and not on the detailed 
abundance ratio in the absence of oscillations. Hence, 
we can utilize the results of [15] for the dependence of 
the 7 Li/ 11 B ratio as a function of sin 2 26 i i3, as shown 
in Figure [T] The temperature dependence is given as a 
linear offset of the model calculations for each value of 
sin 2 26 1 i3. All other model parameter dependences can 
be taken as a correction factor to the 15 M Q calculation. 
That is we write: 

Zi(e 13 ,T u ,R 3a T 12Ca ,M) = [Z(6 13 ) + AZ{T V )\ (10) 

x fma x flil2Ca X /m , 



6 



where Z(9 13 ) + AZ(T V ) is the 16.2 M model calculation 
of [IS]. The temperature correction AZ(T V ) is the top- 
hat distribution of results shown on Figure [T] based upon 
models with different neutrino temperatures |15J. 

The relative ratios as a function of the mixing angle 
should not change much as the progenitor mass changes 
because this will be fixed by the same resonance condi- 
tion independently of the other parts of the star. Hence, 
correction factors for the dependence of the 7 Li/ 11 B ra- 
tio as a function of progenitor mass can be obtained from 
the ^-process (without oscillations) calculations of [T§] . 
These factors are plotted in Figure [2] Similarly, correc- 
tion factors as a function of the R^ a and the R\2Ca as a 
function of progenitor mass can be extracted from |17j . 
These are shown in Figures [3] and [4] Unfortunately re- 
sults in Ref. [17] are only given for progenitor masses 
of 15 and 25 M . For our purposes we have assumed 
constant production factors outside that mass range and 
linearly interpolated production factors inside that mass 
range. 




M 



FIG. 2: Correction factor for the z/-process 7 Li/ 11 B ratio 
relative to the 15 Mq model as a function of progenitor mass 
(deduced from |19|). 



Based upon these assumptions we have preformed 
a numerical evaluation of the five dimensional integral 
given in Eq. ^ . When completed we found that a prior 
equal probability (50/50) for an inverted vs. a normal 
mass hierarchy becomes a 74/26 probability based upon 
the adopted parameters, their prior probabilities, and in- 
fluence on the model results. Hence, although presently 
there are large uncertainties in this approach, based upon 
the present analysis there appears to be a slight prefer- 
ence for an inverted neutrino mass hierarchy. 




FIG. 3: Correction factor for the 7 Li/ 11 B ratio relative to 
the adopted rate (Rz a = 1) as a function of the 3a — > 12 C 
nuclear reaction rate (deduced from [17]). 
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FIG. 4: Correction factor for the 7 Li/ 11 B ratio relative to the 
adopted rate (Ri2Ca = 1.2) as a function of the 12 C(a, 7) ie O 
nuclear reaction rate (deduced from [17]). 



V. CONCLUSION 

In summary, we have shown that an analysis of SiC X 
grains enriched in ^-process material has the potential to 
solve the neutrino mass hierarchy problem in spite of the 
large model uncertainties and the difficulties in determin- 
ing the abundance of 7 Li y and n B„. The fact that 7 Li is 
produced more efficiently than n B in the case of a normal 
hierarchy due to the higher v e energy after oscillations for 
finite 0i3 hints that a combination of accelerator limits, 
SiC grain analysis, and an analysis of supernova models 
might point toward a solution of the neutrino hierarchy 
problem. 
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There are of course, a large number of caveats here. 
For one, we have used a rather crude approximation of 
the model dependence of the 7 Li/ 11 B ratio on various pa- 
rameters from published results in the literature. Clearly 
a more exhaustive self consistent set of model calcula- 
tions would be desirable. Nevertheless, the simple anal- 
ysis applied here at least demonstrates that it might be 
worthwhile to under take such an extensive study. 

Another caveat is that of course there are many un- 
certainties in the meteoritic analysis. For example the 
present analysis assumes that the 7 Li„ and n B„ con- 
densed into the SiC grain in a way that preserves the 
7 Li/ 11 B ratio in the SiC condensation site of the SN 
ejecta. However, because of elemental fractionation be- 
tween the SN ejecta and the SiC grain, the 7 Li/ 11 B ratio 
in the X grain might be different from the true value of 
the ejecta. 

We also note that the analysis here is based upon an 
average of only 7 single X grains and provides only an 
upper limit on the 7 Li 1/ / 11 B 1/ ratio because there is no de- 
tection of (/-process 7 Li (although some individual grains 
do show enhanced 7 Li/ 6 Li [T5]). Both mass hierarchies 
are thus constrained. Nevertheless, we have shown that 



without prior preference for either hierarchy, the inverted 
hierarchy is preferred over the normal hierarchy by a fac- 
tor of ~ 3/1. Clearly, more analysis of ^-process material 
in SiC meteoritic grains is warranted. An analysis of even 
a few more grains could provide a detection of i^-process 
7 Li and substantially improve these limits. Such a de- 
tection would provide both upper and lower limits to the 
7 Li l/ / 11 B l/ ratio and confirm this means to determine the 
neutrino mass hierarchy. 
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